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Abstract 

We explore the properties of test-particle orbits in "bumpy" spacetimes — stationary, reflection- 
symmetric, asymptotically flat solutions of Einstein equations that have a non-Kerr (anomalous) 
higher-order multipole-moment structure but can be tuned arbitrarily close to the Kerr metric. 
Future detectors should observe gravitational waves generated during inspirals of compact objects 
into supermassive central bodies. If the central body deviates from the Kerr metric, this will 
manifest itself in the emitted waves. Here, we explore some of the features of orbits in non-Kerr 
spacetimes that might lead to observable signatures. As a basis for this analysis, we use a family of 
exact solutions proposed by Manko & Novikov which deviate from the Kerr metric in the quadrupole 
and higher moments, but we also compare our results to other work in the literature. We examine 
isolating integrals of the orbits and find that the majority of geodesic orbits have an approximate 
fourth constant of the motion (in addition to the energy, angular momentum and rest mass) and 
the resulting orbits are tri-periodic to high precision. We also find that this fourth integral can 
be lost for certain orbits in some oblately deformed Manko-Novikov spacetimes, leading to ergodic 
motion. However, compact objects will probably not end up on these chaotic orbits in nature. We 
compute the location of the innermost stable circular orbit (ISCO) and find that the behavior of an 
orbit in the approach to the ISCO can be qualitatively different depending on whether the location 
of the ISCO is determined by the onset of an instability in the radial or vertical direction. Finally, 
we compute periapsis and orbital-plane precessions for nearly circular and nearly equatorial orbits 
in both the strong and weak field, and discuss weak-field precessions for eccentric equatorial orbits. 
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I. INTRODUCTION 



The space-based gravitational-wave (GW) detector LISA is expected to detect gravita- 
tional waves generated during the inspirals of stellar-mass compact objects (white dwarfs, 
neutron stars or black holes) into supermassive bodies in the centers of galaxies — extreme- 
mass-ratio inspirals (EMRIs). LISA could detect gravitational waves from these systems for 
several years prior to the plunge of the compact object into the central body and hence ob- 
serve several hundred thousand waveform cycles. Such observations will provide an exquisite 
probe of the strong gravity region close to supermassive central bodies (see [l| for a review). 
In principle, the emitted gravitational waveform encodes the multipole structure of the 
spacetime outside the central object 0]. One of the hopes for LISA EMRI observations is to 
extract this spacetime structure from the data and use it to test whether the central objects 
are indeed Kerr black holes, as we suppose, or something else {2!, 0]. (Intermediate-mass- 
ratio inspirals detectable by Advanced LIGO may reveal the spacetime structure outside 
intermediate-mass central bodies with more modest precision j^.) 

For a Kerr black hole, the spacetime is uniquely determined by the mass and angular 
momentum of the hole and all higher multipole moments depend on these in a simple way 

Ml + iSi = M{ixM)K (1) 

Here Mi and Si are the /'th mass and current multipole moments of the gravitational field, 
M is the mass of the black hole and x is its dimensionless spin parameter, x = Si/M^ = 
a/M. As a consequence of relation ([1]), if the quadrupole or higher multipole moments of 
a supermassive body are measured from an EMRI observation and these are inconsistent 
with the values predicted by its mass and spin, the body cannot be a Kerr black hole with a 
vacuum exterior. The "no-hair" theorem states that, in pure gravity, any pseudo-stationary, 
vacuum and asymptotically flat spacetime containing an event horizon and with no closed 
timelike curves exterior to the horizon must be described by the Kerr metric 0, [6|. If the 
Cosmic Censorship Conjecture is correct, all astrophysical singularities will be enclosed by a 
horizon. It is therefore most likely that the supermassive central bodies which are observed 
to inhabit the nuclei of most galaxies are indeed Kerr black holes. However, LISA should 
be able to test this assumption. Alternatives to Kerr black holes include "dirty" Kerr black 
holes with external masses (e.g., an accretion disk), exotic supermassive stars such as boson 
stars [3], and naked singularities. "Hairy" black hole solutions are also allowed when gravity 
is coupled to other fields, e.g., a Yang- Mills field (these solutions have been shown to be 
uustaMe to perturbations S) 

or a Skyrme field [9|] (stability to generic perturbations is an 
open question). Sufficiently accurate measurements may allow us to distinguish between 
these possibilities. 

In order to prepare us to interpret LISA observations of EMRIs, to identify any deviations 
from Kerr that are manifest in the waveforms and even to facilitate detection of inspirals 
into highly non-Kerr spacetimes, we need to understand how these deviations influence the 
emitted gravitational waveforms. In an extreme-mass-ratio inspiral, the time-scale for the 
orbital inspiral due to radiation of energy and angular momentum is generally much longer 
than the orbital time-scale. We can therefore approximate the inspiral as quasi-stationary. 



2 



by assuming the inspiraling object is always nearly on a geodesic orbit of the spacetime, and 
evolving the parameters determining this geodesic slowly over the inspiral (this is usually 
referred to as the "adiabatic approximation" in the literature 10|] since the fluxes of energy 
and angular momentum used to evolve the sequence of geodesies are computed by assuming 
the object is on an exact geodesic of the spacetime). In this slow-inspiral limit, the emitted 
waveforms depend sensitively on the properties of the geodesic orbits in the spacetime — the 
dominant frequency components in the gravitational waveform at any moment are harmonics 
of the orbital frequencies of the underlying geodesic. We can thus understand some of 
the main consequences of deviations from the Kerr metric by examining the effect of such 
deviations on test particle orbits in the spacetime. By considering a spacetime with an 
arbitrary set of multipole moments, Ryan demonstrated that, for nearly circular and nearly 
equatorial orbits, the periapsis and orbital-plane precessions encoded all of the multipole 
moments at different orders in a weak field expansion 

A multipole moment decomposition is not very practical, however, since an infinite num- 
ber of multipoles are required to characterize the Kerr spacetime. For this reason, Collins 



& Hughes and Glampedakis & Babak [1^ took a different approach and explored test 
particle dynamics in "bumpy" spacetimes, which were constructed as first-order perturba- 
tions of the Schwarzschild and Kerr spacetimes respectively and therefore could be made 
arbitrarily close to Schwarzschild/Kerr by dialing a parameter to zero. Collins & Hughes 
coined the phrase "bumpy" black hole to describe these spacetimes. In their case, the pres- 
ence of stresses exterior to the black hole meant that the horizon could be preserved in the 
presence of the black hole deformation without violating the no-hair theorem. In the present 
case, this name is not strictly applicable since the spacetimes we consider are not black holes 
at all, but rather naked singularities not enclosed by an event horizon. However, the term 
"bumpy" black hole is still a good one to describe how the spacetime appears to an observer 
away from the central object. 

One drawback of the perturbative approach is that the perturbation is not necessarily 
small close to the central body, and so the first-order perturbation theory used to construct 
the spacetime breaks down. As a result, the perturbative solutions may only be used rela- 
tively far from the central object. In this work, we therefore take an alternative approach 
and consider the properties of orbits and inspirals in a family of spacetimes that are ex- 
act solutions of the vacuum field equations of relativity and which include the Kerr and 
Schwarzschild spacetimes in a certain limit. We use a family of spacetimes that were de- 
rived by Manko & Novikov [3] • As exact solutions, the spacetimes are valid everywhere and 
can thus be used to probe the orbital dynamics in the strong-field as well as the weak-field. 
The family has an infinite number of free parameters, which can be chosen to make the mul- 
tipole moments of the spacetime match those of the Kerr spacetime up to a certain order, 
and then deviate at higher order. In this paper, we choose to make the multipole moments 
deviate at the mass quadrupole order and higher, by varying a single parameter, although 
the formalism generalizes to other types of deviation. We use this family of spacetimes as a 
test-bed for an exploration of various observable consequences of deviations from the Kerr 
metric, but we compare to previous work in the literature as we proceed. 

The main new results of the current work are as follows. By studying the properties of 
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orbits in the strong field of tlie spacetime, we find that most geodesies in the spacetime 
appear to have a fourth isolating integral of the motion, in addition to the energy, angular 
momentum and rest mass that are guaranteed by the stationarity and axisymmetry of the 
metric. The corresponding orbits are triperiodic to high accuracy. This was not guaran- 
teed, since the separability of the geodesic equations in Kerr and corresponding existence 
of a fourth integral (the Carter constant) was unusual. Additionally, we find that for some 
oblate perturbations of the Kerr spacetime, there are regions of the spacetime in which there 
appears to be no fourth integral, leading to ergodic motion. If observed, ergodicity would 
be a clear 'smoking-gun' for a deviation from Kerr. Ergodic motion has been found in other 
exact relativistic spacetimes by other authors, although these investigations were not carried 
out in the context of their observable consequences for EMRI detections. Sota, Suzuki and 
Maeda IJ] described chaotic motion in the Zipoy-Voorhees-Weyl and Curzon spacetimes; 
Letelier & Viera 15|] found chaotic motion around a Schwarzschild black hole perturbed by 
gravitational waves; Gueron & Letelier observed chaotic motion in a black hole spacetime 
with a dipolar halo 16|] and in prolate Erez- Rosen bumpy spacetimes 17|; and Dubeibe, 
Pachon, and Sanabria-Gomez found that some oblate spacetimes which are deformed gen- 
eralizations of the Tomimatsu-Sato spacetime could also exhibit chaotic motion 18|]. The 



new features of our current results are the presence of potentially ergodic regions for a wider 
range of magnitudes of the perturbation, and an examination of whether the ergodic regions 
are astrophysically relevant. We find that, in the context of an EMRI, the ergodic regions 
exist only very close to the central body and these regions are probably not astrophysically 
accessible, at least in the Manko-Novikov spacetime family. 

We also look at the properties of the last stable orbit for circular equatorial inspirals. The 
frequency of this orbit will be a gravitational-wave observable, and depends significantly on 
the magnitude of any deviations from Kerr. For certain choices of the quadrupole pertur- 
bation, we find that the last stable orbit is defined by the onset of a vertical instability, 
rather than the radial instability which characterizes the last stable orbit in Kerr. This is a 
qualitative observable that could be another 'smoking-gun' for a deviation from Kerr. 

Finally, we look at the periapsis and orbital-plane precession frequencies. We do this 
primarily for nearly circular and nearly equatorial orbits, since these can be characterized in 
a gauge invariant way in terms of the orbital frequency measured by an observer at infinity. 
Although such precessions were computed by Ryan his results only apply in the weak- 
field. We find results that are consistent with Ryan's in the weak-field, but also explore 
the properties of precessions in the strong-field and find they de pen d significantly on the 



nature and location of the last stable orbit. Collins & Hughes [ll| and Glampedakis & 
Babak [12| did explore strong-field precessions, but they did so as a function of spacetime 
coordinates, rather than as a function of observable quantities which we do here. The 
perturbative spacetimes are also not totally applicable in the vicinity of the last stable 
orbit, so our results are more generally applicable. We also briefly discuss precessions for 
eccentric equatorial orbits in the weak-field and how this is relevant for LISA observations. 

The paper is organized as follows. In Sec. [Til we introduce our chosen family of spacetimes, 
describe some properties of these solutions and discuss our approach to computing geodesies 
in the spacetimes. In Sec. IIIII we analyze geodesies in these bumpy spacetimes and use 
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Poincare maps to identify the presence of an effective fourth integral of the motion. We 
show that most orbits are regular and triperiodic, but also demonstrate the onset of ergodic 
motion in certain oblately deformed spacetimes. In Sec. [IV] we find the last stable orbit for 
circular equatorial orbits and discuss its properties. In Sec. |V] we report our results on the 
periapsis precession and orbital-plane precession in these spacetimes. Finally, in Sec. I VII we 
summarize our results and discuss further extensions to this work. This paper also includes 
two appendices, in which we present results demonstrating ergodic motion in Newtonian 
gravity (Appendix [A]) and an expansion of the precessions in the weak-field (Appendix [Bjl . 
Throughout this paper we will use units such that c = G = 1. 



II. BUMPY BLACK HOLE SPACETIMES 



In this section, we briefly summarize the Manko-Novikov metric [13|. This is the test met- 
ric for which we will explore the dynamics of orbits in Sections IIIIllVl The Manko-Novikov 
metric is an exact stationary, axisymmetric solution of the vacuum Einstein equations that 
allows for deviations away from the Kerr spacetime by a suitable choice of parameters char- 
acterizing the higher-order multipole moments. The presence of these deviations destroys 
the horizon, so this is no longer a black-hole spacetime. However, its geometry is very 
similar to that of a Kerr black hole with additional anomalous multipole moments until 
close to the expected horizon location. We choose a subclass of the Manko-Novikov metric, 
parametrized by a parameter f5. For (3 = 0, the metric corresponds to the usual Kerr metric. 



(In the notation of [l3|, our parametrization corresponds to setting a2 = P and = for 
all n 7^ 2). 

This subclass of the Manko-Novikov metric can be described by a Weyl-Papapetrou line 
element in prolate spheroidal coordinates as (cf. Eq. (1) of 13|): 



ds' = -f{dt-u;d<pr + ef-'e'\x'-y') ( ^ + ) +e f-\x' -l){l~y')d4>', (2) 

— 1 1 — 



5 



where (cf. Eqs. (9, 10, 12, 13 of 0): 

/ = e^^A/B, (3a) 

u = 2ke-'^^CA~^ - Akail - a^Y\ (3b) 

e^^ = exp {2-1'^ A{x^ - lY\l - a^)'^, (3c) 

A = (x^ - 1)(1 + ahf - (1 - y^){h - af, (3d) 

B = [x + l + {x- l)abf + [(1 + y)a + (1 - y)bf , (3e) 
C = (x^ - 1)(1 + ab)[b-a- y{a + b)] + (1 - y'^){b - a)[l + ab + x{l - ab)], (3f) 

^ = PR-^P,, (3g) 

7'^^ln^ + ^(P3P3-AP.) (3h) 

, ^>^/^ a:-y+(-l)^-^(x + y) 
+ ^Z^ - 



i=0 ^ 
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a{x,y) = -«exp(^-2/3(^-l + XJ^^;^^jj, (3i) 

Kx, y) = a exp (^2/3 (^1 + ^^^^^^fe^] ] , (3j) 

R ^ {x' + y'-l)'/\ (3k) 
P„ ^ P„(xy/i?) where P„(x) = ^ ( — j {x' - 1)". (31) 

Here k, a, and /? are free parameters which determine the multipole moments of this space- 
time. The first few multipole moments have the following values (we correct a typo in 
Eq. (14) of dij following [loj): 

Mo = k{l + a^)/{l-a^) Sq = 

Ml = Si= -2ae{l + a^)/{l - a^f 

M2 = -k\(5 + 4a2(i + _ a^Y^\ ^2 = 

M3 = ^3 = 4aA;4[^ ^ 2a2(i ^ ^2)(i _ a^Y^\l{\ - a^). 

(4) 

Therefore, for a given choice of mass M = Mq, spin x = Si/M"^ and anomalous (addi- 
tional to Kerr) dimensionless quadrupole moment q = — (M2 — M^*^") / , the three metric 
parameters are: 



« = ^ ' ^ = ^TT^' ^ = ^ir- (5) 

A given choice of M, x and q uniquely defines the metric. With this definition of q, a choice 
q > represents an oblate perturbation of the Kerr metric, while g < represents a prolate 
perturbation. A spacetime is oblate if it has M2 < 0, e.g., for Kerr M2 = —y^M^. When we 
say a prolate/oblate perturbation we mean a perturbation that makes the spacetime more 
prolate/oblate relative to Kerr. In particular, for —x^ < q < the spacetime is still oblate, 
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although it has a prolate perturbation relative to the Kerr metric. We note that taking 
g 7^ changes all higher moments from their Kerr values, so these solutions deviate not 
only in the mass quadrupole moment but also in the current octupole moment, the mass 
hexadecapole moment etc. 

To present our results, we find it useful to display them in terms of cylindrical coordinates 



p, z and 0. These are related to the prolate spheroidal coordinates x, y by (191 

p = A;(x'- 1)^/2(1 -2/2)^/2^ z = kxy, (6) 
and the line element in cylindrical coordinates is 

ds" = -f{dt - ujd(j)y + [e^\dz^ + dp") + p2rf</.2] . (7) 



A. Spacetime Properties 

The Manko-Novikov spacetimes are vacuum and have the multipolar structure given in 
Eq. (jl]). As a consequence of the no-hair theorem, the spacetimes must therefore either 
lack an event horizon or contain closed timelike curves exterior to a horizon. In fact, both 
of these statements are true. The central singularity is enclosed by a partial horizon at 
coordinates p = 0, \z\ < k. However, this horizon is broken in the equatorial plane by a 
circular line singularity atx = l, ?/ = 0(p = 2: = 0) j2o|. For x = the spacetime is 
otherwise regular, but for x 7^ 0, the spacetimes contain both an ergosphere and a region 
where closed timelike curves exist. The structure of the spacetimes is quite similar to that of 



the (5 = 2 Tomimatsu-Sato spacetime, as described in [21[. The boundary of the ergosphere 
is determined by the condition gu = 0. Inside this region, timelike observers cannot be at 
rest. Such a region is entirely physical, and also exists in the Kerr spacetime, where it is 
of interest since it allows energy extraction via the Penrose process. We show the location 
of the ergosphere for x = 0.9 and various choices of q in the top panel of Figure [H The 
shape of the ergosphere is more complicated when q ^ 0, having a multiple lobed structure. 
This structure is also qualitatively different depending on the sign of g — for g > there 
are three separate ergoregions, one of which intersects the equatorial plane, one which is 
entirely above the equatorial plane and one which is entirely below; for g < there are only 
two regions, one of which is entirely above the equatorial plane and one of which is entirely 
below. 

For a metric of this type, the region where closed timelike curves (CTCs) exist is deter- 
mined by the condition g^^ < 0. In the bottom panel of Figure [H we show the portion of the 
spacetime where CTCs exist for the same choices of g and x = 0.9. Particles orbiting inside 
the CTC region are moving backward in time. This is not inconsistent with relativity, but 
CTC zones are sometimes regarded as unphysical. A spacetime with no CTC zone can be 
constructed by adding an inner boundary in the spacetime, and just using the portion of 
the Manko-Novikov solution exterior to that boundary. 

The CTC zone again has a multiple lobed structure and is different depending on the sign 
of g. We note in particular that for g < the ergosphere does not intersect the equatorial 
plane, although the CTC region does. For g > both regions intersect the equatorial plane, 
and the outermost edge of the CTC region is inside the ergoregion. 



7 




FIG. 1: Spacetime structure for % = 0.9. The upper row shows zeros of gu for q = —1 (left column), 
q = Q (middle column) and q = 1 (right column). This defines the boundary of the ergoregion of 
the spacetime. The region with gu > is shaded. The bottom row shows points where gfj,^ changes 
sign for the same values of q, and the region where g^^ < is shaded. This defines the region 
where closed timelike curves exist. The middle bottom panel is empty since there is no such region 
in the Kerr spacetime. The shape of the two boundaries is qualitatively the same for other values 
of q with the same sign, although both regions grow as \q\ is increased. 



B. Geodesic Motion 



Geodesic motion in an arbitrary spacetime is described by the second order equations 
where the connection coefficients F^^ are given by 

= ^g"^ + - 9p^,f.) ■ (9) 

The spacetimes we are interested in are axisymmetric and time-independent and the metric 
correspondingly has two ignorable coordinates — t and 0. There are therefore two constants 
of geodesic motion: the energy E and the ^-component of angular momentum L^, which are 
given by 

E = -Qui - gt4><i>, = gt^i + g^^cf), (10) 

where a dot ' denotes the derivative with respect to proper time r. Another first integral of 
the motion can be obtained from conservation of the rest mass of the orbiting particle: 

- 1 = gapx'-x^. (11) 
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In practice, we numerically integrate the second-order geodesic equations ([8j) rather than 
use these first integrals, and we use the constancy of E, and Qapx'^x^ as cross-checks 
to verify the quality of our numerical results. The results reported below typically show 
the conservation of these quantities to a few parts in 10^° over the time of integration; see 
Fig. [21 We compute the connection coefficients analytically from expressions for the metric 
functions /, uj and 7 defined in Eqs. (E]). The only difficulty arises at points where a metric 
component g^j^^ vanishes and its inverse g^^ diverges. When this occurs, we analytically factor 
out the terms that tend to zero to avoid issues in numerical integration. To perform the 
numerical integration we write the coupled system of four second-order ordinary differential 
equations ([H]) in first-order form and integrate numerically in C++ via the Bulirsch-Stoer 
method. 
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FIG. 2: The fractional errors in energy E (solid line), angular momentum (dashed line), and 
the quantity gapx^x^ (dotted line) accumulated over 1700 orbits of a geodesic with E = 0.92 and 
Lz = 2.5M in a spacetime with spin x = 0.9 and anomalous quadrupole moment q = 0.95. 



Some general properties of geodesic motion can be understood by using the first inte- 
grals (|TO i) -(fTT i) . The energy and angular momentum conservation equations ( fTOj) can be 
used to write t and in terms of E, Lz, p and z: 



t 



—Egttf, — Lzgtt 



9t4, ~ 9tt94,<i> gt^ - 9tt9<t><i, 

These expressions can be substituted into Eq. (ITT!) to give 

jr-, ip' + = IT—, - - u{p, z) Ef-l^ V^E, Lz, p, z) 



(12) 



(13) 
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The motion in p and z may thus be thought of as motion in the effective potential Vce- In 
particular, since the left hand side of Eq. (fT3|) is strictly positive or zero, motion can only 
exist in regions where > 0. Finding the zeros of the effective potential therefore allows 
us to find allowed regions of the motion. As an illustration, we show the zeros of the effective 
potential in Figure [3] for the simple case of the Kerr metric with spin parameter x = 0-9) 
energy E = 0.95 and angular momentum = 3M. There are two regions of allowed motion 
— one region at larger radius that corresponds to bound orbits, and another region at very 
small radii that corresponds to rising and plunging orbits. 




2 4 6 8 10 12 14 

p/M 

FIG. 3: Effective potential for geodesic motion around a Kerr black hole, with E = 0.95, = 3M 
and X = 0-9- The curves indicate zeros of the effective potential. Allowed orbits are found in the 
small region around p = 0, z = (rising and plunging orbits) or in the region containing p = 10, 
z = (bound orbits). 

We now turn our attention to the Manko-Novikov spacetime with q ^ 0. For spacetimes 
with X = 0, and for spacetimes with x 7^ and g < (prolate perturbation of the Kerr 
metric at large radii), the addition of the perturbation does not fundamentally change the 
nature of the effective potential - there are still two bounded regions, one attached to the 
origin corresponding to rising and plunging orbits and one at larger radii corresponding to 
bound orbits. The shapes of these regions change as |g| is increased and if |g| is increased 
sufficiently at fixed E and the two regions merge, so that all allowed orbits can reach 
the origin. Even after this has occurred, there appear to be two types of orbit in the 
single allowed region - those that rise and plunge and those that undergo many periods 
of radial oscillation. We don't know if the latter remain non-plunging forever in principle. 
In practice, perturbations due to external material or radiation reaction may cause bound 
orbits to diffuse onto plunging orbits over time. For fixed q < 0, the two allowed regions also 
change shape as the energy and angular momentum are varied. In particular, the plunging 



10 



region connected to the partial horizon at p = 0, |2;| < develops a multi-lobed structure. 
For sufficiently large |g| and sufficiently low E and Lz, two of these lobes can touch in the 
equatorial plane. This leads to the existence of circular, equatorial orbits that are unstable 
to vertical perturbations, which we will encounter again in Section IIVI 

For X 7^ ^"^^ q > (oblate perturbation of the Kerr metric at large radii), the behavior 
is qualitatively different. For any arbitrarily small an additional allowed region appears 
in the effective potential, which is bounded away from p = and therefore corresponds to 
bound orbits. For small |g| this new region is very close to p = 0. The other two allowed 
regions still exist, and merely change shape as the value of \q\ is increased. The additional 
bound region is always outside the region where closed timelike curves (CTCs) exist, and is 
therefore in the portion of the spacetime that can be regarded as physical. However, in the 
plane z = the outermost edge of the CTC region touches the innermost edge of the region 
of bound motion. This additional region also extends inside the spacetime ergosphere. 

We consider as an example the case with x = 0.9 and q = 0.95. The zeros of the 
effective potential Ves are plotted in FigurelUfor geodesies with energy E = 0.95 and angular 
momentum = 3M. In this figure there are three distinct allowed regions as described 
above: (i) a foliated "plunging" region connected to p = 0, where all orbits rapidly plunge 
through the horizon (this region also intersects the CTC region); (ii) an inner bound region, 
which is located between p/M ^ 0.72 and p/M ^ 2.12 for the chosen values of E and 
Lz] and (iii) an outer bound region between p/M ^ 2.39 and p/M ^ 13.6. We show the 
trajectory of a typical orbit in the outer region. This has a regular pattern or intersections 
throughout the (p, z) plane, which is characteristic of an orbit with an approximate fourth 
integral. 

If |g| is increased from the value shown in Figure HI the two regions of bound motion 
eventually merge. When this ffist occurs, the "neck" joining the regions is extremely narrow. 
Geodesies exist which can pass through the neck, but this requires extreme fine tuning. As 
|g| is further increased, the neck gradually widens and eventually disappears. At that stage, 
the single allowed region for bound orbits has a similar shape to the outer region of Figure HI 

These general properties of the effective potential seem to be common to all spacetimes 
with g > and x 7^ 0. More relevant for the EMRI problem is to fix q and x ^i-nd to vary 
E and L^. For E = 1 and sufficiently large Lz, there are two regions of allowed motion 
bounded away from the origin, in addition to the plunging zone connected to p = 0, \z\ < k. 
The outermost of the allowed regions stretches to infinity and contains parabolic orbits. The 
inner region of bounded motion is the analogue of the inner bound region described above 
and lies very close to the central object. If the angular momentum is decreased, while keeping 
E = 1, the two non-plunging regions get closer together and eventually merge to leave one 
allowed region that stretches to infinity. For fixed E < 1 the behavior is qualitatively the 
same, except that for Lz^ M there is no outer region (there is a maximum allowed angular 
momentum for bound orbits of a given energy, as in the Kerr spacetime). As Lz is decreased, 
the outer region for bound motion appears and then eventually merges with the inner region. 
Decreasing Lz further eventually causes the bound region to merge with the plunging region. 
At fixed Lzi if there are two distinct non-plunging allowed regions for E = 1, these regions do 
not merge as E is decreased, but the outer region eventually disappears (there is a minimum 
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FIG. 4: Effective potential for geodesic motion around a bumpy black hole with x = O-Si Q. = 0.95, 
E = 0.95, and = 3M. The thick dotted curves indicate zeros of the effective potential. The 
trajectory of a typical geodesic in the outer region is shown by a thin curve. The regular pattern of 
self-intersections of the geodesic projection onto the p—z plane indicates (nearly) regular dynamics. 

allowed energy for orbits of a given angular momentum, as in the Kerr spacetime). If there 
is only one non-plunging region for E = 1, then as E is decreased, this region eventually 
splits into two allowed regions, and the outer region eventually disappears as E is decreased 
further. The properties are similar for all x 7^ 0, but decreasing x with the other parameters 
fixed tends to bring the two allowed regions of motion closer to merger with one another. 



III. ISOLATING INTEGRALS 

The isolating integrals given by the conservation equations f fTO|) - f|TT]) do not completely 
describe the motion, since the motions in p and z are coupled. Thus, solution of the geodesic 
equations requires use of the second order form of those equations (IHl). However, it was 



demonstrated by Carter [22| that in the Kerr spacetime there is a fourth isolating inte- 
gral for geodesic motion, the Carter constant, which arises as a constant of separability of 
the Hamilton- Jacobi equation and was later shown to be associated with a Killing tensor 
in the spacetime. Carter found the form of all metrics that were both Schrodinger and 
Hamilton- Jacobi separable. Imposing the further requirement that the metric be a solution 
of the vacuum Einstein-Maxwell equations leads to the Kerr metric as the only spacetime 
of this form that does not include a gravomagnetic monopole. Thus, the separability of the 
equations in Kerr is somewhat fortuitous and we would not expect that the fourth integral 
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would be preserved when we add an anomalous quadrupole moment as we do here. As a 
consequence, the properties of geodesies might be expected to be somewhat different, and 
might even be ergodic. As mentioned in the introduction, ergodic geodesic motion has been 
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found in other relativistic spacetimes by several other authors 

A fourth integral of the motion essentially gives another relationship between and 
P. Combining this with the effective potential equation fll3l) allows us to eliminate for 
instance and hence obtain an expression for as a function of p and z only. Similarly we 
can obtain an expression for as a function of p and z. 

A standard way to examine equations of motion and look for ergodicity is to plot a 
Poincare map. This involves integrating the equations of motion and recording the value of 
p and p every time the orbit crosses a plane z = constant. From the preceding arguments, if 
a fourth integral exists, the value of p will be a function only of p and z (the function could 
be multi-valued, depending on the order at which the velocities appear in the constants of 
motion). Therefore such a map must show a closed curve. Similarly, if the Poincare map of 
an orbit shows a closed curve for every value of z, then this defines a relationship between 
p, p and z which is then an effective fourth integral of the motion. The Poincare analysis 
thus provides a means to identify whether an effective fourth integral exists or the motion 
is apparently "chaotic". In the latter case, the absence of the integral would be manifested 
on the Poincare maps as space-filling trajectories rather than closed curves. 

The absence of a full set of isolating integrals does not necessarily mean that all orbits 
will exhibit full-blown chaos. For some initial conditions, orbits may show obvious signs 
of ergodicity, while for other initial conditions in the same spacetime, orbits may appear 
to behave in an integrable fashion, suggesting that an approximate additional invariant 
exists. Although this behavior may appear surprising at first glance, it is consistent with 
the predictions of the KAM theorem and with many known examples of chaotic behavior. 
(The KAM theorem, due to Kolmogorov, Arnold and Moser, states that if the Hamiltonian 
of a system with a full set of integrals of motion is analytically weakly perturbed, then phase- 
space motion in the perturbed system will be confined to the neighborhoods of invariant 
tori in phase space, except when angle-variable frequencies of the unperturbed system are 
nearly commensurate, in which case motion will be chaotic [23|.) 

As an illustration, we show in Figure [5] the Poincare map for geodesic motion along orbits 
with three different initial conditions in the Kerr spacetime with the same E, and x 
Figure [3l The Poincare maps are all closed curves, consistent with the existence of the fourth 
isolating integral, the Carter constant. In Appendix [A] we present results for motion under 
gravity in a Newtonian quadrupole-octupole potential and demonstrate the existence of both 
regular and ergodic orbits. This example serves to put the relativistic results described here 
in a Newtonian context. 



1. Poincare Maps for the Manko-Novikov spacetimes 

The regularity properties of geodesies appear to be highly correlated with the nature of 
the effective potential as described in the previous section. For spacetimes with x = 
those with x 7^ but g < 0, all orbits appear to be regular, i.e., they show closed Poincare 
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FIG. 5: Poincare map showing dp/dr vs p for crossings of the z = plane for a sequence of orbits 
in the outer ahowed region of the Kerr spacetime with E = 0.95, = 3M and x = 0-9. The 
closed curves indicates the presence of a fourth isolating integral, which we know to be the Carter 
constant. 

maps similar to those in Figure O These are the spacetimes in the Manko-Novikov family 
that have effective potentials which are qualitatively the same as the Kerr case. 

For q > 0, the effective potential can have two allowed regions for bound motion. What is 
striking is that, whereas orbits in the outer allowed bound region (which corresponds to the 
allowed region in the q = limit) appear to be regular, with closed Poincare maps, those in 
the inner allowed region appear chaotic. In Figures [6] and [7] we show Poincare maps for one 
orbit in each of the outer and inner regions of the effective potential illustrated in Figure H] 
[q = 0.95, E = 0.95, = 3M, x = 0.9). Orbits in the outer region show closed Poincare 
maps, suggesting that the motion is regular or very nearly so and has an approximate fourth 
invariant of the motion. This is reinforced by the projection of the orbit onto the p-z plane, 
which was shown in Fig. |H The geodesic shows a regular grid pattern, with four possible 
velocities at each point, corresponding to ±|p| and ±|i;|. If these orbits do not have a true 
invariant, the regularity of the Poincare map suggests that it may still be possible to find 
an algebraic expression for an approximate constant of the motion. 

Orbits in the inner region, by contrast, seem to fill up all possible points in a subdomain 
of the allowed parameter space (with Ves > 0) and are therefore apparently ergodic in this 
subdomain. It seems hkely, in view of the KAM theorem, that all orbits in the spacetime are 
strictly speaking chaotic, and no true isolating integral exists, but in the outer region there 
is a quantity that is nearly invariant along the orbits 0] . Either the thickness of the region 
mapped out by the chaotic motion is small, or the time over which ergodicity manifests itself 
is very long. From an observational standpoint, whether the motion is actually regular or 



14 



whether only an approximate invariant exists is irrelevant, since the time-scale over which 
ergodicity would manifest itself in the waveform would be much longer than the time during 
which the orbiting object moves on an approximate geodesic. 

It is unusual, given that chaotic and nearly regular regions are generally interspersed in 
most KAM theorem applications 23|, that we find the family of geodesies is divided into 
two distinct regions such that geodesies in one region are ergodic while those in the other 
exhibit nearly regular orbital dynamics. We have been unable to find any strongly ergodic 
geodesies in the outer region, or any non-ergodic geodesies in the inner region. As described 
in the previous section, adjusting the orbital parameters can cause the two allowed regions 
to merge. When this first occurs, the two regions are connected by a very narrow neck. 
The narrowness of the neck means that extreme fine tuning is required to get a geodesic to 
pass through the neck. By choosing initial conditions in the neck, and integrating forwards 
and backwards in time, we obtained orbits that traversed the neck once and found that 
the motion was apparently ergodic while in the inner region, but apparently regular in the 
outer region. This behavior is consistent with the predictions of the KAM theorem, but 
observationally the fact that the orbits in the outer region are technically ergodic does not 
matter as long as they appear regular on long time-scales. We were unable to find an orbit 
that traversed the neck more than once. Further adjustment of the orbital parameters causes 
the neck to widen and eventually disappear. At that stage, most of the orbits appear to be 
regular, but orbits that pass very close to the inner edge of the merged region (i.e., close to 
the CTC zone) have not been fully investigated. 



An alternative explanation of these results [2J] is that the geodesic equations are nu- 
merically unstable in the inner region, and therefore small numerical round-off errors in the 
integration routines are driving the orbits away from their true values. Once again, this 
distinction is not relevant observationally. An astrophysical system harboring an EMRI will 
not be isolated. The gravitational perturbations from distant stars etc. will serve the same 
role in perturbing the orbits as numerical errors might on a computer. The end result — 
that the orbit is apparently ergodic — is the same. 



2. Frequency Component Analysis 

The above conclusions are supported by a frequency- domain analysis of the p and z 
motion in the two regions. The absolute values of Fourier transforms of p{t) and z{t) are 
plotted in Figures M and [9l Fig. [9] shows an absence of clearly identifiable frequency peaks 
for geodesies in the inner region, a result consistent with full-blown chaos. By contrast, 
Fig. [8] shows discrete frequency peaks in the outer region. Generally such frequency peaks, 
corresponding to harmonics of a few fundamental frequencies, occur in problems with a full 
set of isolating integrals. We find that the frequency components measured for the p and z 
motion in the outer region can be represented as low order harmonics of two fundamental 
frequencies at a high level of precision (1 part in 10^ for the first ~ 10 harmonics). This 
multi-periodicity of the geodesies implies that the gravitational waveforms will also be multi- 
periodic. Indeed, we find that an approximate gravitational waveform, constructed using 
a semi-relativistic approximation for the gravitational-wave emission (as used to construct 
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FIG. 6: Poincare map for a geodesic in the outer region of Fig. [H 



Kerr EMRI waveforms in 25|]), is also tri-periodic (the third frequency arises from the (f) 
motion since the observer is at a fixed sky location). The absolute value of the Fourier 
transform of the component of this gravitational waveform is also plotted in Fig. [8] 

and is clearly multi-periodic. This periodicity has important consequences for data analysis 
and parameter extraction. 



3. Comparison to Other Results 



Our results are consistent with previous work by other authors who have found chaotic 
geodesic motion in various spacetimes. Generally, chaotic motion only occurs in the strong- 
field region close to the central object, and for a limited range of geodesic parameters. As 
an example, Gueron and Letelier 17| found chaos in a prolate Erez-Rosen spacetime, which 
represented a deformation of a Schwarzschild black hole. They demonstrated that, for a 
particular value of the energy and angular momentum, when the deformation parameter had 
a value /c2 = —5, there was a single allowed region of bounded motion, but for k2 = —5.02 
the region split into two separate regions. After the split, orbits in the inner region appeared 
chaotic while those in the outer region appeared regular. For the merged region, orbits that 
passed into the inner part also appeared ergodic while those that were purely in the outer 
part looked regular. This is qualitatively very similar to what we have found in the Manko- 
Novikov spacetime, although we find chaotic motion only when x 7^ 0, while Gueron and 
Letelier presented examples for both a perturbed non-spinning black hole and a spinning 
black hole. As a test of our codes, we repeated Gueron and Letelier's calculation and found 
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FIG. 7: Poincare map for a geodesic in the inner region of Fig. HI 

consistent results. As well as providing another example of chaos for relativistic geodesies, 
the results here show some new features. In particular, the inner allowed region appears for 
any g > and as far as we have been able to ascertain the motion is always ergodic in that 
region. This contrasts to the spacetime considered by Gueron and Letelier, in which chaotic 
motion exists only for a small range of k2 (by the time k2 has increased to k2 = —5.1, the 
motion is no longer apparently ergodic). Previous authors have also not considered the issue 
of accessibility of the ergodic region to stars, and we discuss that in the next sub-section. 



Sota et al. discussed what might cause chaos in relativistic geodesic motion, and 
suggested that it might arise either due to a change in the signs of the eigenvalues of the 
Weyl tensor, which would lead to "local instability" or due to the presence of homoclinic 



orbits. The Manko-Novikov spacetimes do contain homoclinic orbits, but Sota et al. [14 
found that this only led to chaos in non-reflection symmetric spacetimes, so this explanation 
probably does not apply here. We have not explored the properties of the eigenspace of the 
Weyl tensor for these spacetimes, but "local instability", could be a plausible explanation 
for our results. The CTC region of the Manko-Novikov spacetime might also be causing 
the ergodicity. The region where ergodic motion occurs touches the CTC region at a single 
point, so the singular behavior of the metric as the CTC region is approached might explain 
the observed behavior, either by causing a region of "local instability" or through some other 
mechanism. 

We note that in the regime where chaos occurs, the perturbation to the Kerr metric cannot 
be regarded as purely quadrupolar, but the deviations in the higher multipole moments are 
also significant. This is similar to the Newtonian result described in Appendix El since 
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FIG. 8: Absolute values of the Fourier transforms of p{t) (solid line), z{t) (dashed line), and the 
gravitational wave component hj^{t) (dotted line) in the frequency domain for an orbit in the outer 
region of Fig. HI The frequency is displayed in units of 1/M; the amplitude scaling is arbitrary. 



we find chaos in the Newtonian quadrupole-octupole potential but not a pure quadrupole 
potential. The relativistic results are somewhat different, however, since we find chaos only 
for X 7^ 0, so for these spacetimes we also need a non-zero current dipole moment to observe 
chaotic behavior. 



4- Accessibility of the Ergodic Domain 

While the existence of ergodic motion is mathematically interesting, an important ques- 
tion for EMRIs that has not been addressed so far is whether ergodicity could ever be 
observed in nature. In other words, is it possible, during the course of an inspiral, for a 
captured object to find itself on an ergodic geodesic? 

In typical astrophysical scenarios, the inspiraling compact object will start out far away 
from the central body with energy close to 1 Unless the angular momentum is very 
small (which in the Kerr spacetime would represent an object on a plunging orbit), this 
will correspond to an orbit in the outer region of allowed motion if two regions exist, so the 
orbit will initially be regular. As the star inspirals, the energy and angular momentum will 
gradually change and this causes the separation between the outermost point of the inner 
region of bound motion and the innermost point of the outer region, Ap, to change. For 
example, when E = 0.99 and Lz = 4.33M in a Manko-Novikov spacetime with x = 0.9, and 
q = 0.95, we find that Ap/M ^ 6.4. When E = 0.95 and = 3M in the same spacetime, 
the separation between regions is only Ap ^ 0.27M. For sufficiently small choices of energy 
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FIG. 9: Absolute values of the Fourier transforms of p{t) (solid line), z{t) (dashed line), and the 
gravitational wave component ^+(t) (dotted line) in the frequency domain for an orbit in the inner 
region of Fig. HI The frequency is displayed in units of 1/M; the amplitude scaling is arbitrary. 



and angular momentum (e.g., E = 0.92 and = 2.5M) only a single region remains. This 
suggests that the two regions will come closer together as energy and angular momentum are 
radiated away during an inspiral, until they eventually merge. We conjecture that d{Ap)/dt 
is always negative; that is, the two regions are always merging rather than separating. To 
test this conjecture, we must explore the behavior of Ap along an extreme-mass-ratio inspiral 
characterized by slowly evolving E and L^. 

To do this, we use an approximate scheme to evolve the energy and angular momentum 
during an inspiral. Our scheme is based on combining exact relativistic expressions for the 
evolution of orbital elements with approximate post-Newtonian formulae for energy and 
angular-momentum fluxes. This scheme was previously devised to describe EMRIs into 
Kerr black holes 26|] and has been shown to give reliable results in that context. For the 
current calculation, we must augment the fluxes with an additional post-Newtonian term 
to represent the effect of the anomalous quadrupole moment q on the evolution of energy 
and angular momentum. A Kerr black hole has quadrupole moment M2/M^ = — x^- It is 
the quadrupole moment that leads to the lowest order terms in in the expressions for 
the energy and angular momentum radiated during an inspiral. Therefore, to include the 
excess quadrupole moment, we just change the terms in the flux expressions to -|- q, 
while leaving the lower order terms unchanged (this approach was also used in 27|). We 
then numerically find the roots of the effective potential Vcs = in the equatorial plane at 
various times and compute the evolution of Ap along the inspiral. 

The result of one such computation of Ap is plotted in Fig. [TOl That figure corresponds 
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to an inspiral in a spacetime with x = 0-9) ^"^^ Q = 0.95. The inspiral starts out at 
p = lOOM with an orbital inchnation of 60 degrees and initial eccentricity e = 0.8 (these 
orbital parameters correspond to ~ 0.9982 and ~ 5.0852M) and proceeds until plunge. 
The separation between the inner and outer bounded regions gradually shrinks, until the 
two regions merge (on the plot, this is shown as Ap = 0). Afterward, the bounded regions 
remain joined until eventually merging with the plunging region. 




0.5 1 1.5 2 2.5 3 3.5 



t/M ^^06 

FIG. 10: The evolution of the separation Ap between the inner and outer bounded regions in 
the equatorial plane along an inspiral in a Manko-Novikov spacetime with x = 0.9 and q = 0.95. 
Ap = means that the two regions have merged and there is a single bounded region. 

We have found the same qualitative behavior described above for a wide range of param- 
eter choices. Therefore, in all these cases, our conjecture is true — the inspiraling object can 
never find itself in the isolated inner region where all orbits appear to be ergodic. We should 
point out, however, that we have carried out this numerical investigation only for a range 
of specific choices of and initial orbital parameters, and have used an approximation 
to the energy and angular momentum radiated during an inspiral. This is therefore not a 
definitive proof that chaotic motion can never be observed in the course of an inspiral in the 
Manko-Novikov spacetime. 

Assuming this evolution really is typical, there are two important consequences. Firstly, 
an inspiraling object can never end up in the inner of two allowed regions of bound motion, 
where ergodic motion is prevalent. Secondly, inspirals always start out in a phase where the 
motion is regular. This is very important, since it will allow the systems to be detected in 
this early inspiral stage by gravitational-wave detectors using matched filtering or a time- 
frequency analysis. The inspiraling object will eventually end up in the merged region formed 
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after the two regions of bounded motion converge. Both ergodic and regular geodesies exist 
in that region, so in principle the particle could find itself on an ergodic orbit. However, 
most orbits in the merged region appear to be regular so it would require fine tuning to put 
the object onto such a geodesic (e.g., the "neck traversing" geodesies discussed earlier). It 
thus seems unlikely that this would occur in practice. 

Although these results apply only to the Manko-Novikov family of spacetimes, the con- 
clusions are consistent with other examples of chaotic geodesies in relativity. For instance, 
in the prolate Erez- Rosen spacetime considered in [13], if an object had arrived in the region 
where ergodic motion is observed during the course of an inspiral, its orbital energy and 
angular momentum would have been larger earlier in the inspiral. However, if either the 
energy or angular momentum is increased from the values that give ergodic motion, the 
effective potential changes so that it has only one allowed region, which includes "escape 
zones" connected to the central singularity. All geodesies in such a zone plunge into the 
central object in a short time so an astrophysical inspiral could not persist through that 
zone. We deduce that for that spacetime as well the ergodic region is inaccessible to objects 
captured at large distances. 

If there was some other mechanism that could put an inspiraling object onto an ergodic 
geodesic, there is the question of how the ergodicity could be identified in practice. Detection 
of EMRIs will rely on matched filtering or possibly time- frequency techniques In either 
case, it will probably not be possible to identify the gravitational radiation as being emitted 
from an ergodic orbit, but only that radiation from a regular orbit has ceased. It is clear from 
Figure [9] that during an ergodic phase, the emitted power is spread among many harmonics, 
which will consequently not be individually resolvable. This radiation will increase the 
broadband power in our detector, whereas if the orbit had plunged the radiated power would 
rapidly die away. However, the energy released during a typical EMRI is comparatively 
low, so it is unlikely that we could identify the presence of such broadband power over 
the instrumental noise. Therefore, the chances are that we will not be able to distinguish 
observationally between an inspiral that "ends" at a transition into an ergodic phase and 
one which ends by plunging into a black hole. 

One potentially observable signature of ergodicity would be an inspiral that turned "off" 
and "on" as it progressed through ergodic phases interspersed with regular phases. This 
would occur if the object could move into and out of the inner ergodic region during an 
inspiral, but the preceding analysis indicates that this shouldn't happen. An object on a 
"neck-traversing" geodesic would also show this behavior. However, the periods where the 
orbit is ergodic serve to randomize the phase of the orbit in the regular periods. A signal 
of this type would only be observable if each apparently regular phase could be individually 
resolved with enough signal-to-noise ratio. This would require a very narrow "neck" in order 
to trap the orbit for many cycles in the regular zone. However, fine tuning of the energy 
and angular momentum is necessary to make the neck very narrow, so if an object was on 
such an orbit, the neck would be widening rapidly as energy and angular momentum were 
radiated away. In practice, it is doubtful that sufficient signal-to-noise would accumulate to 
allow a detection to be made before the neck widened too much. 

We conclude that, for astrophysically relevant inspirals in the Manko-Novikov spacetime 
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family, an object would probably not end up on an ergodic geodesic. If some other mech- 
anism conspired to put an object on such an orbit, it is unlikely that we would be able 
to identify this in gravitational-wave observations. If these findings carry over to a more 
generic class of spacetimes, then chaotic motion is merely a mathematical curiosity which 
is unlikely to manifest itself practically or be important for gravitational-wave data analysis 
considerations. 



IV. LAST STABLE ORBIT 



During an inspiral into a Kerr black hole, an EMRI will evolve quasi-stationarily through 
a sequence of near-geodesic orbits as orbital energy and angular momentum are radiated 
away. There is a minimum energy (which is dependent on angular momentum) for which 
bound orbits exist. When the inspiral reaches that separatrix, the object will rapidly plunge 
into the central body. The gravitational radiation emission undergoes a transition at this 
point, and so the frequency of this last stable orbit is in principle another quantity that is 
observable from the detected gravitational waves. For a Kerr inspiral, the 'transition' is a 
rapid die-off in the gravitational-wave emission as the particle plunges into the black hole. 
If the central object is not a black-hole, the radiation may persist for longer after the last 
stable orbit is passed [3], but there will still be a significant qualitative change in the emitted 
radiation as the orbit changes suddenly at that point. We focus on the innermost stable 
circular equatorial orbit in this analysis, since this is well defined in these spacetimes. 



A. Circular equatorial orbits 

The geodesic equations for an arbitrary spacetime ([8]) may be written in the alternative 
form 

d / dx"\ 1^ dx^dx" 

= T^d^gua (14) 



dr V dr J 2^"''" dr dr 

For a circular, equatorial orbit in an axi- and reflection-symmetric spacetime of the form 
d?]), dp/dr = dz/dr = d^p/dr^ = 0; hence the p-component of the geodesic equation (fn|) 
gives 

dpg^xpcj)^ + 2dpgt^i(t) + dpgut^ = 0, (15) 

in which a dot denotes d/dr as before. We can thus express the azimuthal frequency as 
observed at infinity = (p/t in the form 

= 9;:^, • 

where the +/— signs are for prograde and retrograde orbits respectively. In the equatorial 
plane, the right-hand side is a function of the spacetime parameters and p only, so given a 
particular choice of azimuthal frequency Q^p, Eq. (fT6l) can be inverted to determine the value 
of p such that a circular orbit at that p has frequency Q^. 
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Equation (ITTll provides another relation between t and 4>, from which we can deduce 

_ 1 

i = {-9tt - 2n^gt^ - VL^g^^) ^ , (17) 

and then the energy and angular momentum equations (ITOl) give us E and Lz as a function 
of p for circular equatorial orbits. 

B. Innermost Stable Circular Orbit 

The location of the innermost stable circular orbit (ISCO) in the equatorial plane can 
be found using the effective potential (fT3l) . Circular equatorial orbits are located at simul- 
taneous zeros and turning points of Ves, where Vcs = dVcs/dp = dVca/dz = 0. As we will 
see in Section |V] the second derivatives of V^s determine the frequencies of small oscillations 
about the circular orbit. For the circular orbit to be stable, we need the orbit to sit at 
a local maximum of Kg, i.e., we require d'^Vcs/dp'^ and d'^Vcs/dz'^ to be negative. In the 
following we will use V^p(p) (V^^(p)) to denote the value of d'^V^s/dp'^ {d'^Ves/dz'^) evaluated 
for the circular equatorial orbit at radius p. For the Kerr spacetime, Kz(p) < at all radii, 
but Vpp{p) has a single root at a critical radius pisco- This tells us that the orbit becomes 
radially unstable at that point, which defines the ISCO. For x = 0, pisco ~ 4.90M , while for 
X = 0.9, Pisco ~ 1.25M for prograde orbits and pisco ~ 7.705M for retrograde orbits. Note 
that p is a cylindrical Weyl coordinate, which is why these results differ from the familiar 
black-hole ISCO radii, which are normally quoted in Boyer-Lindquist coordinates. 

For the Manko-Novikov solutions with x = 0, the shape of the functions Vpp(p) and 
Vzz{p) does not change significantly as q is increased with g > 0: Vzz{p) < everywhere 
and Vpp{p) = has a single solution that defines the ISCO. However, as |g| is increased 
with q < 0, there is a transition in behavior at g ~ —0.163. For q < —0.163, the function 
Vpp{p) has two zero-crossings. Thus, in addition to the radially-stable circular orbits at 
large radii, we find additional such orbits exist very close to the central singularity. If 
is increased still further, the two roots converge at g ^ —0.654 and for q < —0.654 radially 
stable orbits exist at all values of p. However, at the point where the second branch of 
the radial roots appears, there is also a transition in the shape of Vzz{p), so that there 
are now orbits which are vertically unstable. For q < —0.163, the ISCO is defined by this 
vertical instability, rather than the radial instability characteristic of the Kerr spacetime, 
and Manko-Novikov spacetimes with g > 0. In the range —0.654 < g < —0.163, there are 
two regimes where stable circular orbits exist — an outer zone with p > pisco? and an inner 
zone with pisco < P < Posco ("we use "OSCO" to indicate "outermost stable circular orbit" 
and Pisco to denote the ISCO for the inner set of circular orbits). The energy and angular 
momentum of an orbit at the "OSCO" are greater than the energy and angular momentum 
at the ISCO of the outer zone, pisco- Thus, an object inspiraling from large distances on a 
circular equatorial orbit will reach pisco and plunge into the central body, rather than finding 
itself in the inner range of circular orbits. Compact objects could only find themselves in 
the inner range if they came in on an eccentric/inclined orbit and then radiated away energy 
and angular momentum in exactly the right proportions. It is therefore unlikely that this 



23 



o 

u 



Radially unstable - first branch 
Radially unstable - second branch 
Vertically unstable 



0.22 
0.2 
0.18 
0.16 

O 

M 0.14 



s 



0.12 
0.1 
0.08 
0.06 
0.04 





Radially unstable - first branch 




Radially unstable - second branch 




/ Veitically unstable 

; 


- 


; 


/ 

/■ 


■■/ 
•■1 

I 
/ 

- 







q q 
FIG. 11: Properties of the equatorial ISCO in spacetimes with x = 0, as a function of q. We show 
the p coordinate of the ISCO (left panel) and the dimensionless frequency of the orbit at the ISCO 
(right panel). As described in the text, the ISCO radius has three branches, depending on whether 
it is determined by one of the two branches of radial instability or the branch of vertical instability. 
These branches are indicated separately in the diagram. For values of q where all three branches 
are present, the dashed line denotes the "OSCO" and the dotted line denotes pisco as discussed 
in the text. Allowed orbits lie above the curve in the left panel, and below the curve in the right 
panel. 



inner zone would be populated in practice. However, any object on a circular equatorial 
orbit in this inner zone would reach pisco and then plunge into the central body. 

In Figure [11] we show the location of the ISCO as a function of q for spacetimes with 
X = 0. We also show the orbital frequency at the ISCO as a function of q, computed using 
Eq. (fT6|) . For spacetimes with spin, the behavior is qualitatively similar, but there are now 
two ISCO radii, corresponding to prograde and retrograde orbits respectively. We show 
results for a spin of x = 0.9 in Figure [T2l We note that the ISCO radius is always outside 
the boundary of the causality-violating region of the spacetime. For x 7^ and g > 0, 
the ISCO radius is determined by the energy at which the outer allowed region for bound 
motion (which is a single point for a circular equatorial orbit) merges with the inner allowed 
region. In that case when the object reached the ISCO it would undergo a transition onto 
an eccentric/inclined geodesic. 

The value of the ISCO frequency depends not only on q but also M and x- However, 
as we shall discuss in the next section, it is possible to measure these other parameters 
using precessions measured when the orbit is in the weak-field. Thus, the ISCO frequency 
is a powerful probe of the nature of the spacetime since it can be very different even for 
comparatively small deviations from Kerr. 
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FIG. 12: As Figure [TTl but now for a spin of x = 0.9. There are now two ISCO curves, one for 
prograde orbits and one for retrograde orbits. The ahowed range of orbital frequencies is given by 
the region in between the two curves in the right hand panel. 



V. PERIAPSIS AND ORBITAL-PLANE PRECESSIONS 



In Section IIIII we saw that astrophysically relevant orbits in the Manko-Novikov space- 
time are multi-periodic to high precision. In such cases, there is no smoking-gun signature 
that indicates the presence of "bumpiness" in the spacetime. Instead, the imprint of the 
spacetime bumpiness will be observationally apparent in the location of the last stable orbit, 
as discussed in the previous section, and in the following ways: (1) in the three fundamen- 
tal frequencies of the gravitational waves generated while the inspiraling object is on an 
instantaneous geodesic orbit; (2) in the harmonic structure of the gravitational-wave emis- 
sion, i.e., the relative amplitudes and phases of the various harmonics of the fundamental 
frequencies; and (3) in the evolution of these frequencies and amplitudes with time as the 
object inspirals. A full analysis of the accuracies that could be achieved in observations 
would involve computing gravitational waveforms in the bumpy spacetimes, performing a 
Fisher-Matrix analysis to account for parameter correlations, and comparing to a similar 
analysis for Kerr. That is beyond the scope of this paper. However, we can examine the 
first of these observational consequences by comparing the fundamental frequencies between 
the bumpy and Kerr spacetimes. 

The complication in such an analysis is to identify orbits between different spacetimes. 
Identifying orbits by the p and z coordinates is not gauge- invariant, since the meaning of 
these coordinates depends on the spacetime structure. Identifying orbits via the energy 
and angular momentum is gauge-invariant, but these quantities are not directly measurable 
observationally. However, circular orbits in the equatorial plane of the spacetime are char- 
acterized by a single observable — the azimuthal frequency of the orbit. We can use this 
frequency to identify circular equatorial orbits in different spacetimes. 

Precession frequencies are absent in exactly circular equatorial orbits. However, if the 
circular orbit is perturbed radially, it will undergo small oscillations at the radial epicyclic 
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frequency, which is characteristic of the periapsis precession frequency at that radius. Like- 
wise, if the orbit is perturbed vertically it will undergo small oscillations at the vertical 
epicyclic frequency, which is characteristic of the orbital-plane precession frequency at that 
radius. We thus compare these epicyclic frequencies, as a function of the circular orbital 
frequency, between Kerr and bumpy spacetimes. This comparison was employed by Ryan, 
who used it to derive his theorem stating that all spacetime multipole moments are encoded 
in the gravitational waves generated by nearly-circular, nearly-equatorial EMRIs j^. 

An eccentric equatorial orbit can be characterized by two observables — the orbital 
frequency and the periapsis precession frequency. These two frequencies can therefore be 
used to identify orbits in different spacetimes (provided there is an orbit with corresponding 
frequencies in the Kerr metric). Likewise, the orbital-plane precession frequency can be used 
to identify inclined orbits between spacetimes^. With such an identification, differences in 
the multipole structure of the spacetime will show up only in the relative amplitudes of the 
harmonics and in the evolution of the fundamental frequencies over the inspiral. We will 
discuss this some more at the end of this section, but a full analysis requires treatment of 
inspiral in an arbitrary spacetime and is beyond the scope of the current paper. 

A. Epicyclic frequencies 

The frequency of epicyclic motion can be derived by perturbing a circular, equatorial 
orbit in either the radial or vertical direction. The second order geodesic equations (HM for 
z and p take the form 




Here X denotes either p or z. The dependence on dt/dr and d</)/dr can be eliminated by 
using the energy and angular momentum conservation equations to express these in terms 
of E, Lz, p and z, as in Eq. ( fT2l) . Using this form of the equations we can take a circular, 
equatorial orbit, p = pc, z = 0, and perturb it either in the radial direction, p = pc + Sp, 
z = 0, or in the vertical direction, p = pc, z = 6z. Considering the equations of motion at 
leading order in the small orbital perturbation, it is easy to see that the frequencies of these 



The 'orbital-plane' is not well defined in the strong field. However, we know the gravitational waves 
should be triperiodic and, in the weak-field, the three periods correspond to the orbital period and the 
two precessions. When we refer to 'orbital-plane precession frequency' we really mean the frequency 
component of the orbit that corresponds to orbital-plane precession in the weak-field. This will be the 
frequency of the vertical motion, averaged over many orbits. 
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small epicyclic oscillations are given by 
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(19) 



As before, X denotes either p (for the radial epicyclic frequency Qp) or z (for the vertical 
epicyclic frequency fi^). The same result can be derived starting from the effective potential 
equation (|T3ll : the frequency fix is determined by d'^Ves/dX'^ evaluated at the circular orbit. 

B. Precessions 

We are interested in precessions rather than the epicyclic frequency. We define the periap- 
sis precession as the number of cycles by which the periapsis advances per radial period (i.e., 
over one complete epicyclic oscillation). Likewise, the orbital-plane precession is defined as 
the number of cycles by which the azimuthal angle to the highest point of the orbit advances 
during one vertical oscillation. These precessions, which we denote by px, are related to the 
epicyclic frequencies, fix, by 

Px = ^- 1- (20) 

\ix 

The behavior of the precessions can be understood in terms of what happens in the weak- 
field, far from the black hole, and in the strong-field, close to the ISCO. In the weak-field it is 
possible to derive expressions for the precessions as functions of the orbital frequency. This 
was originally done for nearly circular, nearly equatorial orbits by Ryan who demon- 
strated that the various spacetime multipole moments enter the precession rate expansion 
at different orders of (Mfl^)'^. This was the basis for a theorem that, in principle, the weak- 
field precessions can be used to extract the lowest order spacetime multipole moments. The 
weak-field expansion of the precessions is summarized in Appendix [Bl 

In the strong-field, we find that one or the other precession diverges at a certain frequency. 
This frequency corresponds to the frequency of the ISCO. To understand what is happening, 
we use the effective potential (fT3|) and consider radial oscillations. For the energy and angular 
momentum corresponding to the circular equatorial orbit at radius p = pc, the effective 
potential in the equatorial plane takes the form Veff(/0, z = 0) = —V {p){p — p^){p — pcY . Here 
V{p) is a function that is strictly positive for p > p^. The radius p_ is the other solution to 
Kff(p, z = 0) = 0, and p_ < pc- As the ISCO is approached, the effective potential develops 
a point of inflection at the location of the turning point rather than a maximum since 
P- Pc- The epicyclic frequency for radial oscillations is fi^ oc V{pc){pc — P-), which thus 
tends to zero as the ISCO is approached. The corresponding periapsis precession diverges. 
The radius p_ corresponds to an unstable circular orbit, and associated with any unstable 
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circular orbit is a bound, eccentric orbit that has an infinite period — the object comes 
in from larger radii, and asymptotically approaches the radius of the circular orbit. This 
is referred to as a "homoclinic" orbit, or as a "zoom-whirl" orbit in the EMRI literature. 
As the ISCO is approached, a small perturbation from the location of the circular orbit 
will put the object onto an orbit that is close to the homoclinic orbit associated with the 
unstable circular orbit. Hence, it takes a very long time for the object to complete a radial 
oscillation, but it is moving rapidly in the azimuthal direction the whole time, building up 
a large periapsis precession. 

This understanding leads us to expect the precession to diverge at the location of the 
ISCO, and this divergence should be like (pc — Pisco)"^^^, or (fi(/,,isco — f^^)"^'^^- The above 
argument applies to an ISCO defined by a radial instability (as in the Kerr metric). As 
we saw in Section llVt the ISCO in the Manko-Novikov spacetimes can be determined by 
the onset of a vertical instability. In that case, the above argument still applies, but it is 
now the orbital-plane precession that will diverge as the ISCO is approached. This provides 
another potential 'smoking-gun' for a deviation from the Kerr metric. The divergence in 
the precession at the ISCO arises as a result of one of the two epicyclic frequencies going to 
zero. It is these frequencies that will in principle be observable in the gravitational waves. 
If an inspiral is observed starting in the weak-field and up until the last stable orbit (LSO), 
the different frequency components could be tracked, and one frequency will tend to zero 
as the LSO is approached. This is in principle an observable, and if it is the orbital-plane 
precession that goes to zero the central body cannot be a Kerr black hole. A more careful 
treatment of the gravitational-wave emission will be required to understand how practical 
it will be to make such observations. 

In Figures [TSHM] we show the precessions as a function of MQ^ for a variety of values of 
q. In Figures [T7l - [T9l we present the same results, but now we show the differences between 
precessions in a bumpy spacetime with a given q and precessions in the Kerr spacetime with 
the same spin parameter x'- ^Px = Px — p^'^'^'^- The variable Apx represents the number 
of cycles of difference, so for instance a value of App = 0.1 means that the orbits in the 
two spacetimes, although having the same azimuthal frequency, would drift an entire cycle 
out of phase in the epicyclic radial oscillation within ten radial orbits. We do not show 
results for the difference in the orbital-plane precession for x = 0? since there is no orbital- 
plane precession in the Schwarzschild spacetime and hence that plot would be identical to 
Figure fT5l 

Figures [13] and [H] show the periapsis precession rp(f2^) for x = while Figures [TH 
and [18] show the periapsis precession for x = 0.9. We see that as the value of q de- 
creases from zero, the periapsis precession decreases relative to the corresponding value 
in the Kerr /Schwarzschild spacetime. By contrast, if q is increased from zero, the periapsis 
precession increases. In spacetimes with non-zero spin, the difference is more extreme for 
prograde orbits than for retrograde orbits. This is presumably because retrograde orbits do 
not get as close to the central object, and so do not "feel" the strong-field deviations in the 
bumpy metric. 

For q > —0.5, the radial epicyclic frequency np(fi^) approaches zero as the ISCO is 
approached and the periapsis precession rp goes to infinity for the reasons described above. 
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FIG. 13: Periapsis precession pp versus azimuthal frequency $7^ for x = and various values of q. 



This is not true of the q < —0.5 spacetimes shown, since for those the ISCO is defined by 
a vertical instabihty. Figure [TS] shows the orbital-plane precession r^i^^) for x = and 
Figures [16] and [19] show the orbital-plane precession for x = 0-9- As for the case of the 
periapsis precession, the orbital-plane precession behaves qualitatively differently depending 
on the sign of q. The orbital-plane precession is greater for g < and smaller for g > 
compared to the non-bumpy value. As expected, the orbital-plane precession tends to a 
constant at the ISCO for q > —0.5, while it diverges for q < —0.5, since the ISCO for the 
latter spacetimes is defined by a vertical instability as discussed earlier. 

Previous authors have looked at precessions in "bumpy" spacetimes. As mentioned above, 
Ryan 0] derived a weak- field expansion for the precessions. Collins & Hughes 11 1 looked 
at precessions for eccentric equatorial orbits in a perturbed Schwarzschild spacetime, and 
Glampedakis & Babak 12| did the same for a perturbed Kerr black hole. However, both 
pairs of authors did this by comparing orbits with the same coordinates, which is rather 
unphysical. Our results are consistent with this previous work in the weak-field, but our 
calculation is the first that can be applied in the strong field, since Ryan's work used a 
weak-field expansion, and the other work used perturbative spacetimes that break down 
close to the central body. The behavior in the approach to the ISCO is thus a new result. 

It is possible to fit the precessions as a sum of a weak-field expansion (as given in Ap- 
pendix 



plus a term A/ a/^0,isco — ^4,- However, only a comparatively few weak-field 
terms are required to give a good fit, implying that the divergence at the ISCO limits the 
number of multipole moments that can be recovered from such an expansion. To quantify 
this statement properly, we must do a full analysis, that includes the effect of inspiral, uses 
an instrumental noise curve to restrict the observable bandwidth and accounts for param- 
eter correlations via the Fisher Matrix. We can do this by constructing semi-relativistic 
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FIG. 14: Pcriapsis precession pp versus azimuthal frequency $7^ for x = 0.9 and various values of 
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inspiral waveforms for bumpy spacetimes in the same way that has been used for Kerr in- 



spirals 25|, l26[. This is beyond the scope of the present paper. However, there are several 
things that we can take away from the current results — the location of the ISCO has a 
strong influence on precessions that could be observable, in particular the nature of the in- 
stabihty that deflnes the ISCO could be a clear indicator of a non-Kerr system; precessions 
can be very different in the strong fleld in the presence of a deviation; circular orbits with 
frequencies very different from the Kerr value exist in some bumpy spacetimes, so another 
observable signature would be that an inspiral persists at frequencies inside the Kerr ISCO. 



C. Effect of eccentricity 

As discussed above, the measurement of the precessions as a function of orbital frequency 
for nearly circular, nearly equatorial orbits would in principle allow measurement of the 
spacetime multipole moments 0] . In practice, however, the precessions will only be manifest 
in the observed gravitational waves if the orbit is not circular and equatorial, so we need 
to understand how the dependence of the precessions on azimuthal frequency differs when 
we relax the assumption of near-circularity. In the following, we shall focus on the periapsis 
precession of eccentric but equatorial orbits. 

The eccentricity of an orbit modifles two things — 1) the frequency associated with the 
periapsis precession as a function of the orbital frequency; 2) the relative amplitudes of 
different harmonics of these two frequencies in the observed GWs. To accurately compute 
the dependence of the harmonic structure on eccentricity for a generic orbit, we need to 
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FIG. 15: Orbital-plane precession pz versus azimuthal frequency for x = various values of 
q. We do not show the case q = here, since there is no orbital-plane precession in Schwarzschild. 



know details of GW generation in a spacetime with arbitrary multipole moments. This 
is a difficult problem, so we focus on the effect of eccentricity on the periapsis precession 
frequency itself. We consider an eccentric equatorial orbit in the Kerr spacetime, and use 
flfP to denote the average azimuthal frequency (i.e., the total advance in over one radial 
period, divided by the period of the radial motion). We define an orbital eccentricity, e, 
such that the ratio of the Boyer-Lindquist radii of the periapsis, Vp, and apapsis, r^, of the 



radial motion is rp/r^ 



e)/{l + e). With these definitions, the periapsis precession as 



defined above can be found to be 



Pp 



-4x 



-(34-18e2)x 



e2) 2 



+ 



3(18-7e2 + 2x 
4 




l-e2 



(21) 



where we are expanding in the weak field, MQ^ <^ 1. The corresponding result for a 
spacetime with an excess quadrupole moment can be found at lowest order by replacing the 
term in with + Q, since the quadrupole moment of a Kerr black hole is discussed 
earlier. 

In the circular limit, e = 0, the expansion ( !2TI) allows us to extract M from the coefficient 



of ^'^^^ , X from the coefficient of fi^, q from the coefficient of Q'^'^ etc. However, if we 
expand to lowest order in the eccentricity, e, it is clear that the effect of a small excess 
oblate quadrupole moment g > could be mimicked, at leading order, by an eccentricity 
evolving as = 2{Mil^)~'^^^q. The two possibilities are then distinguished by knowing how 
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FIG. 16: Orbital-plane precession pz versus azimuthal frequency for x = 0.9 and various values 
of q. 



the eccentricity should evolve with MQ^. 

The expansion (12T]) contains redundancy, since the coefficient of MQ^J^ also depends only 
on the lowest current moment, x- If the eccentricity of the orbit did not evolve with time the 
ffist four terms in the expansion would determine M, x, q and the eccentricity e, and higher 
terms would determine the remaining multipole moments as in the circular case. However, 
the eccentricity does evolve with time. In practice, we will only observe EMRIs as they 
evolve through a finite range of frequencies (determined by the detector sensitivity). During 
that period, the evolution will be driven entirely by gravitational-wave emission. This means 
that we can quantify the eccentricity of the orbit by a single number — the periapsis at which 
the eccentricity was equal to 1 if we integrated the inspiral backwards in time, assuming 
a purely GW driven inspiral. Specifying this parameter and the multipole structure of the 
spacetime determines the eccentricity as a function of MQ^. Determining this relationship, 
however, requires knowing the details of GW emission in an arbitrary spacetime. 

A complication arises because the ratio Mf2^/(1 — e^)^/^ tends to a constant at the point 
where e = 1. Assuming that this occurred in the weak-field, Mf2^/(1 — e^)^/^ <^ 1, this 
can be seen by considering the leading order term in de/dQ^ in the weak-field (see for 



instance 



3) 



de 



-(304 + 121e2)(l -e2)e 



d(Mfi^) 3{Mn^) (96 + 292e2 + 376^) ' 
Denoting X = 1 — e^ and expanding in the limit MQrp — 0, X — > 1, we find 
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FIG. 17: Difference between periapsis precessions in a bumpy spacetime with x = ^^'^ the 
Schwarzschild spacetime, App(Q0, q) = pp{^l^, q) — pp{^l(i,, q = 0). 



in which Xq is a constant that is related to the periapse at "capture" when X = 0. If the 
capture occurs in the strong field, the ratio Mil^/{1 — e^)^^'^ would still tend to a constant if 
we integrated backward until e ^ 1. Although the inspiral would not be observed as e — > 1, 
that section of the inspiral does affect the portion that we can observe. 

We now substitute the asymptotic behavior (j23l) into Eq. ( j2Tll . to obtain an expansion 
of the periapsis precession as a function of the angular frequency in the form Pp = ao + 
a2{Mil^Y^^ + as{Mil^) H — ■ , where Oq, a2, 03 etc. are constants. In contrast to the circular 
case, each of these coefficients depends on all the spacetime multipole moments, so multipole 
extraction from the periapsis precession expansion is no longer straightforward. The reason 
for this qualitative difference between circular and eccentric orbits is that it is only possible 
to observe an eccentric inspiral over a finite range of periapsis, since the orbit is captured 
with a certain finite periapsis, while a circular orbit could inspiral from infinity. The various 
multipole moments have different radial dependencies, thus if one can observe the precession 
frequency at any radius it makes sense that all the moments can be separately extracted, 
while this is more difficult if only a finite section of the spacetime is explored. 

In practice, this difficulty also arises when observing a circular inspiral, since the radiation 
can only be detected in a certain frequency range. One can parameterize an observation by 
the frequency at the start of the observation, Qq = ^(i)it = 0) . A Taylor series expansion of 
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FIG. 18: Difference between periapsis precessions in a bumpy spacetime with x = 0.9 and the Kerr 
spacetime with x = 0.9, App(0<^, q) = Pp{Q.^, q) - Pp{^^, g = 0). 



the precession (see Eq. flBlSp in the Appendix) then gives 

Pp = (3 (Mfio)^ - 4x (Mfio) + ^ (9 + + g) (Mfio)^ + ■ ■ 

+ {2 (Mno)^ - 4x (Mfio) + 2 (9 + + g) (Mfio)^ + ^ ~ 



no 
2 



In this kind of expansion the multipole moments again contribute at all orders. However, 
provided the initial frequency MQq <^ 1, the dominant piece of the constant term, bo, is 
{MQq)^ , so this term can be used to estimate M. Similarly, the dominant piece of 2bo — 3bi is 

4x {MQq), so this can be used to estimate and that estimate of x can be used to improve 

4 

the estimate of M from bo. The dominant piece of bo — bi + 862 is (9 + + q) /2 (Mf2o)^, 
so this can be used to estimate the excess quadrupole moment q and so on. In the same 
way, if an eccentric inspiral is observed in a regime where the initial frequency is small (and 
hence the frequency at capture was also small), we can use the same type of expansion and 
use combinations of the coefficients to successively extract each multipole moment and the 
initial eccentricity. To do this requires an expansion of — eg as a function of Q^/Qo — 1. 
The necessary derivatives de"^ /d{MQ^) are known in the weak-field, and to lowest order 



in the multipoles (see, for example, reference 26|). However, this calculation is somewhat 
involved and beyond the scope of this paper. 
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FIG. 19: Difference between orbital-plane precessions in a bumpy spacetime with x = 0.9 and the 
Kerr spacetime with x = 0.9, Apz{i}^, q) = Pz{^}^, q) - Pzi^cp, q = 0). 

The above discussion indicates that the periapsis precession rate can be used on its own to 
measure the multipole moments from an eccentric equatorial inspiral, although this is more 
difficult than for the circular equatorial case. However, the value of the precession is not 
the only observable. As mentioned earlier, the relative amplitude of the various harmonics 
of the orbital frequencies is a powerful probe of the orbital eccentricity. To exploit this 
harmonic structure we also need to know how the amplitudes of the harmonics depend on 
the spacetime multipole moments. However, if the deviations from the Kerr metric multipole 
structure are small, we could imagine using the Kerr harmonic amplitude relation to estimate 
the eccentricity (and inclination) of the orbit, and then use the precessions to extract the 
multipole structure. Proper insight can be gained using the approximate semi-relativistic 
waveforms described earlier or post-Newtonian expansions of the gravitational waveforms. 
Such an investigation will be an important extension of the current work. 

VI. SUMMARY 

In this paper we have discussed various observational signatures that could leave an im- 
print on an EMRI gravitational waveform if the spacetime in which the EMRI was occurring 
deviated from the Kerr metric. We have seen that some orbits in "bumpy" spacetimes lack 
a fourth integral of the motion, and appear ergodic. Geodesies in the Kerr spacetime have a 
complete set of integrals, so if an apparently ergodic orbit was observed it would be a clear 
signature of a non-Kerr central object. However, regions of ergodic motion only appear very 
close to the central object, in a regime which is probably inaccessible to a star inspiraling 
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from large distances. Most astrophysically relevant orbits are regular and appear to possess 
an approximate fourth integral of the motion, and the orbits are tri-periodic to high accu- 
racy. The deviations of the central body from Kerr then manifest themselves only in the 
changes in the three fundamental frequencies of the motion and the relative amplitude of 
the different harmonics of these frequencies present in the gravitational waves. For nearly 
circular, nearly equatorial orbits, the dependence of the precession frequencies on the orbital 
frequency is well fit by a combination of a weak field expansion that encodes the multipole 
moments at different orders, plus a term that diverges as the innermost stable circular orbit 
is approached. The frequency of the ISCO and its nature (whether it is defined by a radial 
or vertical instability) is another observable signature of a non-Kerr central object. 

To derive these results, we have focussed on a particular family of spacetimes due to 
Manko and Novikov [l3|. However, we expect the generic features of the results in the weak 
field and as the ISCO is approached to be true for a wide range of spacetimes. Chaos has been 
found for geodesic motion in several different metrics by various authors Q, H M, E3, 



In all cases, however, the onset of chaos was qualitatively similar to what we found here 
— it occurred only very close to the central object, and for a very limited range of orbital 
parameters. The conclusion that gravitational waves from ergodic EMRIs are unlikely to be 
observed is thus probably quite robust. 

Precessions for spacetimes that deviate from the Kerr metric have also been considered 
by several authors 2], 
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Our results agree with this previous work in the weak-field 
as they should. However, the results in the present paper are the first that are valid in 
the strong-field since previous work was either based on a weak-field expansion j3] or a 
perturbative spacetime 
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The main feature of the precessions in the strong-field — 
the divergence of one of the precessions as the ISCO is approached — is expected from 
spacetime-independent considerations and therefore should be a g eneral feature of inspirals 
in any spacetime. The present work, and earlier research ll|, [l^, has also considered only 
solutions that first differ from the Kerr metric in the mass quadrupole moment. The Manko- 
Novikov solutions 13|] include spacetimes that first differ at higher orders. While we have 
not considered such solutions, we expect the generic features to be similar. The precessions 
will be closer to the Kerr values for a greater fraction of the inspiral, and the ISCO will be 
at a different frequency, but the qualitative behavior in the approach to ISCO should be the 
same. 

The next step in understanding how gravitational- wave detectors might identify non-Kerr 
central objects from EMRI observations is to consider the gravitational waveforms produced 
during an inspiral. Any analysis should account for both parameter correlations and the 
finite bandwidth and observation time of g ravitational-wave detectors by using a Fisher Ma- 
trix analysis. Glampedakis and Babak 12] constructed approximate gravitational waveforms 
generated by orbits in a perturbed Kerr spacetime, but they considered only waveforms from 
geodesies (i.e., not inspirals) and compared waveforms with the same orbital parameters. 
These are not observable quantities (unlike the frequency of the orbit which we used as a ba- 
sis for comparison here) and such a calculation does not account for parameter correlations. 
Barack and Cutler 27|] did a full Fisher Matrix analysis of this problem, and estimated 
that a LISA observation of an EMRI could measure the quadrupole moment of a body to 
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an accuracy of 10~^ while simultaneously measuring the mass and spin to 10~^. That cal- 
culation was based on an approximate waveform model devised to describe Kerr inspirals. 
The expressions governing the inspiral were modified by adding the leading order effect of a 
quadrupole moment to the energy and angular momentum fluxes. The waveform generation 
part of the algorithm was left unchanged. Although this result is a good guide, the calcu- 
lation contained a number of inconsistencies. For Kerr inspirals, semi-relativistic "kludge" 
waveforms based on combining exact geodesic motion with approximate gravitational-wave 



emission formulae have proven to give accurate results 25 



26 . The same method could 



be used to produce waveforms for inspiral in the Manko-Novikov spacetimes, by changing 
the geodesic equations and augmenting the inspiral fluxes appropriately. Such an approach 
will not generate totally accurate gravitational waveforms, but it will reproduce the main 
features of the orbit — the precession frequencies, the orbital shape and the frequency of the 
ISCO. A study of gravitational waveforms generated in "bumpy" spacetimes will provide 
useful guidance for future detectors such as to what precision an observation could deter- 
mine that an inspiral is an inspiral into a Kerr black hole and how well observations can 
distinguish different types of deviation from Kerr, e.g., an exotic central object from a naked 
singularity from a Kerr black hole with external matter. 
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APPENDIX A: CHAOTIC MOTION IN NEWTONIAN GRAVITY 

The classic example from astrophysics of a system that exhibits chaos in classical (New- 
tonian) gravity is the two dimensional Henon-Heiles potential V{r,6) = r^/2 + r"^ sin(3^)/3 



l^see [28| for example). Gueron and Letelier 16|] also found chaos in the Paczyhki-Witta 
potential ($ = M/(r — rs), where rs = 2M is the Schwarzschild radius) with a dipolar 
perturbation. Neither of these spacetimes is reflection symmetric, so for a better analogy to 
the relativistic spacetimes considered in this paper, we examine the Newtonian quadrupole- 
octupole potential 

Here, M, M2 and M4 denote the monopole (mass), quadrupole and octupole multipole 
moments of the potential. Stationarity and axisymmetry ensure that energy E and angular 
momentum = r^d0/dt are conserved as usual, which leads us to the Newtonian analogue 
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of the effective potential equation (fT3|) 

\ ((1)^-^(^)1 ^ ^-(^■^---) 4 - 1) - f - 

where we have replaced the standard Newtonian energy by the relativistic expression [E"^ — 
l)/2 for consistency with f|T3|) . The equation of motion in this potential takes the usual 
form d^r/dt^ = — V$. If we take the multipole moments to have the values M2 = 2M^ and 
M4 = lOM^, and choose the angular momentum to be Lz = 1.7M, we find that for a range 
of values of the energy E, bound orbits occur quite close to the origin. For sufficiently large 
values of E, there is a single allowed region for motion (defined by Vcs > 0). Orbits in that 
regime appear to be regular, and show closed Poincare maps. If the energy is reduced, the 
allowed region eventually splits into two separate regions, one bounded away from r = 0, 
and one connected to r = 0. Orbits in the outermost region after this transition exhibit 
ergodic behavior. In Figure [201 we show four plots. Two of these plots are for an orbit with 
E = 0.82, which exhibit regular behavior. The other two are for E = 0.81 and exhibit 
ergodic behavior. We choose the initial conditions of both orbits to be p = = z and 
p = 3M, with z determined from the assigned energy (1A2I) . The upper panels in the figure 
show the orbit in the (p, z) plane, and the boundary of the allowed region of motion (defined 
by Vcs = 0). The lower panels show Poincare maps for the two orbits. The ergodicity of the 
orbit with E = 0.81 is quite evident from the Poincare map. We also find that this orbit 
fills up the entire allowed range of p and z. By contrast, the regular orbit with E = 0.82 
explores only a narrow torus in space. 

A thorough examination of when ergodicity appears in this potential, as a function of 
energy, angular momentum and the multipole moments M2 and M4 is peripheral to the 
focus of this paper. However, the results presented here provide a Newtonian example to 
which we can compare the relativistic results of Section IIIII 



APPENDIX B: WEAK FIELD PRECESSIONS 
1. Relativistic Precession 

In Boyer-Lindquist coordinates, the energy, angular momentum and rest-mass conserva- 
tion equations ( fTOl) - (fTTl) for geodesic motion in the Kerr metric can be used to derive the 



equation of motion in the form (see for instance 29| ) 




dt J ^ 

{E{r^ + a^) - aLzf - A (r^ + (L, - aE)^ + cos' 6 + cos' 9(1 - E^)) 
2 (E((r2 + a2)2 /A - sin' 6) - 2MaL,r/A)^ 



(Bl) 



38 





0.01 



0.005 



Q. 

T3 



-0.005 



-0.01 




0.2 - 



0.1 



Q. 

T3 



2.7 2.75 



2.85 2.9 2.95 



3.05 



-0.1 - 

-0.2 - 
2 



2.5 



3.5 



p/M p/M 

FIG. 20: Example of onset of chaos in the Newtonian quadrupole-octupole potential ()Aip . All 
plots are for orbits which start with /5=0 = z,p/M = 3 and have specific angular momentum 
Lz = 1.7M. The left hand panels are for energy E = 0.82, while the right hand panels have 
energy E = 0.81. The top two plots show zeros of the effective potential, Ves = 0, as defined by 
equation ()A2p . and the paths followed by the orbits in the (p, z) plane. The bottom two plots are 
Poincare maps for crossings of the z = plane in each case. 

where A = — 2Mr + a^, and a = Mx- The prograde equatorial circular orbit at radius r 
has energy and angular momentum 



E 



1 - 2i;2 + av^/M 
^1 -3v^ + 2av'yM 
1 - 2avyM + a^v^/M^ 

rv 



V^l - 3t;2 + 2av^/M 
where = M/r. The frequency of a prograde circular orbit is given by 



(B2) 
(B3) 



d0 
dt 



-3/2 + aVM 



(B4) 
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The epicyclic frequencies for radial and vertical perturbations of the orbit are given by the 
second derivatives of the right hand side of equation (]B1|) with respect to r and 9 (the right 
hand side of Eq. (IBip is the effective potential for the Kerr spacetime). To obtain the form 
of these frequencies in the weak field, we wish to expand in 1/r. With some manipulation 
and keeping terms up to only, we obtain the expansion 

where we use flp, Qz to denote the radial and vertical epicyclic frequencies to be consis- 
tent with the results earlier in the paper. With further manipulation, expressions for the 
precessions, px, as a function of the orbital frequency, Q^p, may be derived 

Pp = 3 (Mn^f^ - 4x (MQ^) + 1(9 + x') {MQ^f^ - 34x (Mfi^)i 

+^{l35 + 67x'){Mn^f + --- (B7) 

Pz = 2x {M%) - -x^ {M%Y + {Mn^f + ■■■ (B8) 

Results for retrograde orbits may be obtained by the substitutions x ^ ~X-, ~^ ~^4> and 
Lz — > —Lz in the above expressions (NB fi^ < for retrograde orbits, so — is equivalent 

to ir^^i). 



2. Precession due to a Quadrupole Moment 

The precession induced by a quadrupole moment can be derived using the Newtonian 
quadrupole potential 



$ = ---^4 1-3% . (B9) 



Here r = -^/x^ + + is the distance from the origin, z is the vertical coordinate and we 



will use p = ^y x"^ + to denote the cylindrical polar radial coordinate. The radial equation 
of motion in this potential takes the form 



2\dt J 2p2 r 2r3 V r"^ , 

and the energy, angular momentum and orbital frequency of a circular, equatorial orbit with 
radius p are 

Q__M 
4p3 2p 



E 




^<P = \-T + 7;-. (Bll) 
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Differentiating Eq. (IBIOI) twice with respect to p and 2;, we find the epicychc frequencies 
take the form 

Hence we derive the precession frequencies 

p. = \§p (Affl*)* + • ■ ■ ■ (B14) 



The lowest order form of these expressions was also given in Collins and Hughes [11| , although 
they expressed the precession in terms of a radial coordinate, rather than the observable 
VL^. We also use a slightly different definition for the quadrupole moment Q so that it is 
consistent with Q = — x^M^ for the Kerr metric. As we would expect, the leading-order 
terms in these expressions agree with the leading order terms in in the Kerr expressions. 

Combining this result with Eq. ( 1B8I) . we obtain the weak- field precessions for the Manko- 
Novikov solution with spin parameter x and excess quadrupole moment q 

Pp = 3 {M%Y - 4X (Mfi^) + ^ (9 + + g) {Mn^Y - 34x 
+^ (135 + 67x' + 39g) {M%f + ■■■ 

Q 

p, = 2x (Mfi^) - 2 (x' + q) {Mn^Y + (8x' - 3g) {Mn^f + ■■■. (B15) 

In the above, the lowest order term that is omitted is the order at which the excess current 
quadrupole moment would first contribute. This result is also given in Ryan [2I, although 
he quotes an expression for fip/fi^, where VLp is equal to — Vtp. Our result is consistent 
with his once this is taken into account. We note that some of the terms in expression ( 1B15I1 



come from relativistic corrections to the effect of the quadrupole moment. These cannot be 
derived using only the results quoted in this appendix, but are given in Ryan's paper 
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